Quantum spin liquid in a kagome lattice spin-1/2 XY model with four-site exchange 
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We study the ground state phase diagram of a two-dimensional kagome lattice spin-1/2 XY model 
(J) with a four-site ring exchange interaction (K) using quantum Monte Carlo simulations. We find 
that the superfluid phase, existing in the regime of small ring exchange, undergoes a direct transition 
to a 1,2 quantum spin liquid phase at (K/ J) c ~ 22, which is related to the phase proposed by Balents, 
Girvin and Fisher [Phys. Rev. B, 65 224412 (2002)]. The quantum phase transition between the 
superfluid and the spin liquid phase has exponents z and v falling in the 3D XY universality class, 
making it a candidate for an exotic XY* quantum critical point, mediated by the condensation of 
bosonic spinons. 



The coveted quantum spin liquid state, proposed more 
than three decades ago by Fazekas and Anderson [1], re- 
mains surprisingly elusive. Although recent experiments 
[2H1] provide tantalizing evidence of their existence, the 
scarcity of spin liquid states in microscopic models be- 
lies their resistance to theoretical characterization. Even 
a basic unifying notion for which microscopic ingredi- 
ents are required to promote spin liquid states in models 
is muddied. Geometric frustration is the main suspect; 
recent Density-Matrix Renormalization Group measure- 
ments have championed the case for a spin liquid as the 
long-debated groundstate of the kagome lattice Heisen- 
berg antiferromagnet [5]. However, the apparent dis- 
covery of a gapped spin liquid state in quantum Monte 
Carlo (QMC) simulations of the honeycomb lattice Hub- 
bard model at half filling seems to contradict some long- 
held notions for where spin liquids might lie [BJ. The 
extensive theoretical framework [7HTU]. developed over 
decades, continues to undergo refinement [TTI - [T3] moti- 
vated by these and other discoveries from large-scale com- 
puter simulations of a relatively small number of models. 

Part of the difficulty in finding models that harbor spin 
liquids in two-dimensions (2D) and higher is the fact the 
important ingredient of frustration typically leads to the 
infamous sign problem and preclusion of the model from 
being studied by scalable QMC techniques. A significant 
development in this front came in 2002, when Balents 
et. al. [14] proposed a sign-problem free Hamiltonian 
of spins on the kagome lattice which gives rise to a Z2 
spin liquid phase. Away from their exactly soluble point, 
the spin liquid appears to be robust, as several mod- 
els containing XY terms and constrained potentials have 
now shown evidence of this spin liquid ground state |15l - 
117] . In this paper, we demonstrate that a Z2 spin liquid 
phase can also be stabilized in a model with competition 
between purely kinetic terms, namely two-site XY and 
four-site "ring" exchange interactions. This Hamiltonian 
is amenable to large-scale quantum Monte Carlo (QMC) 
studies, which allows us to characterize the spin liquid 
state - and the quantum phase transition into this state 
- in great detail. We find that there exists a single di- 




Figure 1: (color online) (a) kagome lattice and a labeling 
convention for the indices of an operator Pijki- Two primitive 
vectors di, a~2 are shown, (b) Three different orientations of the 
bow-tie plaquettes in the kagome lattice, (c) The two spin 
plaquette configurations between which the K term can act 
upon. The up and down arrows correspond to up and down 
spins, respectively. 

rect quantum phase transition into the spin liquid from 
the superfluid phase, characterized by exponents z and 
v falling in the 3D XY universality class. This makes 
the transition an interesting candidate for an exotic XY* 
quantum critical point [18H20] caused by the condensa- 
tion of bosonic spinons. 

We define the standard J-K model by means of the 
U(l) symmetric spin-1/2 Hamiltonian: 

H = -jJ2Zij-Kj2 p m (!) 

(ij) (ijkl) 

where J is the nearest neighbor coupling and K is the 
four-site ring exchange. The bond and plaquette oper- 
ators are given by B tJ = (S+SJ" + S^S+) = 2(SfS? + 
S?SV) and P ijM = (S+STS+Sr + Sr S +S^S+). The 
bond operator Bij represents the standard nearest- 
neighbor XY exchange interaction, where (ij) denotes a 
pair of nearest neighbor sites on a 2D kagome lattice. The 
plaquette operators Pijki (with labeling convention as in 
Fig. [T^,) make up the purely XY part of the full four-site 
ring exchange interaction, where (ijkl) are sites on the 
corners of a bow-tie plaquette. Similar J-K models have 
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proven a fertile ground in the search for exotic phases 
and phase transitions in other lattices. Recent propos- 
als include Exciton Bose Liquids [21] , d-wave Bose liquids 
(DBL) 22], and other exotic gapless Mott insulators [2"5] . 
The square-lattice XY ring-exchange model was the first 
example investigated as a candidate deconfined quantum 
critical point between a superfluid and Valence-Bond- 
Solid (VBS) phase, after QMC investigations pioneered 
by Sandvik et. al. |Mj elucidated the ground-state phase 
diagram. Interestingly, Buchler et. al. [2 5) have recently 
presented the design of a ring-exchange interaction for 
cold bosonic atoms in 2D optical lattices, making the in- 
vestigation of such models more than just a theoretical 
curiosity. Here, we examine the groundstate behavior of 
the kagomc-lattice J-K model with QMC, restricting the 
sign of K to be positive to avoid the sign problem. In 
the limit of a vanishing ring exchange, e.g. K — > 0, the 
J-K model reduces to the spin S —1/2 XY model, where 
the hopping term J drives the system into a superfluid 
phase at low temperature. The other limit, K — > oo, 
leads to a pure ring-exchange model on the the kagome 
lattice, closely related to a model that has been explic- 
itly shown to stabilize a spin liquid phase by Balents et 
al. [14] . In the presence of a fixed nearest neighbor hop- 
ping, one may expect that the ring exchange term drives 
the system from the superfluid state at a small K to & 
spin liquid state at an extremely large K. Below, we show 
that the superfluid and spin liquid phases are separated 
by a single phase transition at (K/J) c = 21.8 (Fig. |2|. 

In order to study the groundstate phase diagram of 
the kagome J-K model using QMC, a variation of the 
stochastic series expansion (SSE) algorithm [26] is used. 
We consider a finite size system mapped to a torus to 
have periodic boundaries in both directions, with size 
Li = nidi and L2 — 712^2 (where di and 02 are primitive 
vectors and ai = a 2 = 1) as depicted in FigjlJ,. The total 
number of sites in the simulation cell is N s — rii x ri2 x 3 
sites. For simplicity, we use ri\ = 712 = L in our next 
discussion. In order to perform the QMC, several essen- 
tial modifications must be made to the SSE algorithm for 
the square-lattice J-K model using the conventional no- 
tations, i.e. J, C and K-vertices. In the kagome case, the 
number of C and K vertices are the same in the SSE 
QMC scheme (see Ref. [55] Fig. 2). However, unlike 
the square lattice model, the J operators are only as- 
sociated with the pair of two opposite spins highlighted 
with the arrows on the bow-ties plaquette as shown in 
Fig. [l^. In other words, the only valid bonds are ij and 
kl (whereas invalid bonds are il and jk) for this partic- 
ular bow-tie. As the results, the number of J vertices 
are reduced to 16, instead of 32 as in the case of square 
lattice. Technically, this leads to a unique and smaller 
set of directed loop equations, although the algorithm 
implementation in both cases are similar. Moreover, like 
the square-lattice model, simulation efficiency depends 
on the inclusion of special multi-branch cluster updates, 
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Figure 2: (color online) Data collapse of the superfluid den- 
sity p s for P/L = 2/J, (K/.J) C = 21.8, z = 1, and v = 0.6717. 
Inset: Schematic ground state phase diagram of the Hamilto- 
nian Q as a function of the 4-site ring exchange K/J. The 
critical value of ring exchange (K/ J) c = 21.8 separates the 
superfluid phase (SF) from the quantum spin liquid phase 
(QSL). 



which facilitate sampling of C — > K [26] vertices in the 
large- K regime of the model. The square-lattice multi- 
branch cluster update must again be modified in a non- 
trivial way to accommodate the kagome lattice geome- 
try, where the reduced set of unique J vertices reduces 
the number of possible vertex transformations. Similar 
to the square lattice model, a single cluster can be con- 
structed and flipped with probability 1 (a Wolff update) , 
or all clusters can be constructed and flipped indepen- 
dently with probability 1/2 (a Swendsen-Wang update) 
|26j . In the regime of large- if, we have observed that 
the Wolff-type update is more efficient to the Swendsen- 
Wang type. 

In order to characterize the various phases in the 
model, we compute the spin stiffness, as well as the spin 
and plaquette structure factors. The spin stiffness (su- 
perfluid density in the boson representation) is defined 
in terms of the energy response to a twist <j) m the pe- 
riodic boundary of the lattice: ps — jj 9 ■ Start- 
ing at small K where the SSE algorithm is known to 
be ergodic, the superfluid density ps decreases with in- 
creasing ring exchange, eventually vanishing at a critical 
value (K/J) c . At this point, the disappearance of su- 
perfluid density signals an insulating phase; the smooth 
change of ps as a function of K/J suggests a contin- 
uous quantum phase transition. To examine the na- 
ture of this transition further, we apply scaling theory: 
ps = L- z F ps {L x l v (K c - K)/J,/3/L z ) [37], where F ps is 
a universal scaling function, L is the linear system size, 
/3 is the inverse temperature or the imaginary time, z is 
the dynamical critical exponent and v is the correlation 
length exponent. At a fixed ratio f3/L z and sufficiently 
large inverse temperature, a plot of psL z as a function 
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Figure 3: (color online) The spin-spin structure factor 5 s (q) 
as a function of (qx,q y ) for the two states: (a) superfluid 
state with Kj J — 16 and (b) insulating state with K/J = 26. 
(L = 12) 
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Figure 5: (color online) The plaquette-plaquette structure 
factor 5 p (q) as a function of {qx,q y ) for the two states: (a) 
superfluid state with K/J = 16 and (b) insulating state with 
K/J = 26. (L = 12) 




Figure 4: (color online) The real-space image of the spin-spin 
correlation function, {S§Sf). The reference spin (black) at 
site is arbitrary. The color scale is from red (SgSf ) = 0.05 
to blue for (SgSf) = —0.05 with increasing transparency as 
we approach (S§ Sf) = (from either direction). The picture 
shows most spins as transparent, or (SqS^) w for most 

of K/J should intersect at the critical point (K/J) c . Us- 
ing this approach, we have identified the intersection at 
(K/J) c = 21.8 using z = 1 and /3/L = 2/J for different 
system sizes. As illustrated in Fig. [2j one may thereby 
determine the correlation length exponent v, by plotting 
psL z as a function of [(K c — K)/ J\L X I V . All the curves 
collapse into a single curve with v — 0.6717 taken from a 
conventional 3D XY universality class [55]. As discussed 
below, because of the condensation of bosonic spinons at 
the transition, it is predicted that this scaling behavior 
actually belongs to the different universality class called 
XY*, in which an anomalous exponent r\ is much larger 
than the one belonging to the 3D XY universality class, 
Vxy* > Vxy ~ 0.04 [IgH2Dl l2g ] . This prediction could 
be addressed in future QMC studies. 

The adherence of the supcrfluid-insulating transition 
to this universality class lends credence to the hypothe- 
sis that the insulating phase is in fact a Z2 spin liquid, as 
seen in previous studies [TSJ [T7] . Further evidence can be 
obtained by examining the z-component of the spin-spin 
correlation functions, which we average over the prop- 
agated states, i.e. {S z k Sf) = i(^Epi 1<T I(p)^f(P)). 



where n is the number of non-identity operators in the 
SSE operator list, the z-component Sf. l — (l/2)ir| l (with 
o\ l — ±1) [55]. These are used to construct the spin 
structure factor S s (q x ,q y ) = £ ^ e i{r k - ri ).^ S z S z ^ 
where k and I are lattice sites; r; = (xi,yi) is the lat- 
tice coordinate (with lattice spacing chosen as unity) 
and q = (q x ,Qy) is the wave vector. In order to ex- 
amine the characteristics of the insulating phase, we in- 
vestigate the spin structure factor S s (q) as illustrated in 
Fig. [3] The spin structure factor looks similar in both 
the superfluid state and the insulating phase correspond- 
ing to K/J — 16 and K/J = 26, respectively. At any q 
vector, there are no extensively-scaling Bragg peaks (re- 
quired for long-range order) in the insulating state. The 
non-extensive features persisting in large system sizes, 
e.g. L — 12, reflect the existence of short-range correla- 
tions. This evidence rules out a possible charge density 
wave with a broken translation symmetry. It is also use- 
ful to visualize the real space structure of the spin-spin 
correlation function in the insulating state: K/J = 26 is 
shown in Fig. [4] The real-space picture shows a local con- 
figuration that reduces nearest neighbor hopping energy 
(J) while the reference spin (black) is also able to partic- 
ipate in exchange with all four plaquettes it is a part of. 
This configuration is in line with what one would expect 
to find through performing local perturbation theory. 

In Fig. [5] we present the plaquette structure factor 
S P (q x ,q y ) = ^E a , b e 4(r »- r " )q (Pan), where P a is the 
plaquette operator with the plaquette subscript a defined 
in Fig. [TJ:. This characterizes the modulations of the 
plaquette operator expectation value (Pijki)- Similar to 
the spin structure factor, peaks in the plaquette structure 
factor indicate a VBS phase only if they survive in the 
thermodynamic limit. Clearly, the plaquette structure 
factor S p (q) shows no evidence of any extensively scaling 
peaks [31j at any q vector, ruling out the possibility of 
having a VBS for the insulating phase. 

In summary, we have studied the ground state phase 
diagram of the kagome lattice spin 1/2 XY model with 
four-site ring exchange model using SSE-QMC calcula- 
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tions. We find a continuous quantum phase transition 
at (K/J) c « 22 from the superfluid phase to a quan- 
tum spin liquid. Finite-size scaling studies confirm to 
high accuracy that the quantum critical point has a dy- 
namical exponent z — 1 and a correlation length ex- 
ponent v = 0.6717 falling into the 3D XY universality 
class. Above this transition, we observe a large region of 
a gapped spin liquid phase, which dominates the phase 
diagram. Spin and plaquette correlation functions are 
observed to have no long-range order over the entire Bril- 
louin zone. 

This spin liquid phase should be adiabatically con- 
nected to the Z2 spin liquid phase proposed by Balents 
et. ah, which is demonstrated to be stable in other re- 
lated models [T5HT7] . This is supported by the similarity 
in the q-space spin structure factor between our ring- 
exchange system and these models. In addition, the fact 
that the superfluid-insulator transition is observed to be 
in the 3D XY universality class also provides strong evi- 
dence that the featureless insulator phase is in fact a spin 
liquid. In future studies, it would be interesting to inves- 
tigate the finite temperature phase diagram accessible 
with SSE-QMC. In particular, smoking-gun evidence of 
the emergent Z 2 topological order should be possible us- 
ing Renyi entropy measurements ; 30] , however due to the 
significantly stronger critical coupling in the J-K model, 
these studies may require more computational resources 
than the model in Ref. (17] . 

This model provides a novel route to studying ex- 
otic criticality, since the scaling behavior of the quan- 
tum phase transition may actually belong to an exotic 
XY* universality rising from the condensation of bosonic 
spinons. Measurements of the anomalous dimension ex- 
ponent would be sufficient to confirm this novel behavior, 
especially since the numerical value expected to be much 
larger (of order r\ = 1.37 [5D]) than 77 at a conventional 
transition. Just as interestingly, extensions of this model 
with symmetry breaking interactions are expected to give 
rise to a measurable vison-confinement, which should also 
be manifest as an anomalously large (but unique) value 
of i]. Other interesting phases, such as exotic supersolid 
phases, may also be present in straightforward extensions 
of the model, such as with the application of a uniform 
external field. 

Our work has identified a broad new class of XY ring- 
exchange model which supports an extended region of 
quantum spin liquid, amenable to large-scale QMC sim- 
ulations without the sign problem. The rise of the spin 
liquid due to the competition of purely kinetic interac- 
tions, which do not contain conventional geometric frus- 
tration, marks significant progress in the understanding 
of what ingredients are necessary to promote spin liquids 
in realistic models [25] . The abundance of possible exten- 
sions to this model which are likely to show exotic phases 
or phase transitions suggests that large-scale simulation 
work on kagome-lattice J-K models has only yet begun. 



We are grateful to S. Isakov, T. Senthil and C. Xu 
for fruitful discussions. R.G.M. would like to acknowl- 
edge the support and hospitality of Microsoft Station Q. 
This work was made possible by the computing facili- 
ties of SHARCNET. Support was provided by NSERC 
of Canada. 



[1] P. Fazekas and P. Anderson, Philos. Mag. 30, 423 (1974). 

[2] J. R. et ai, Phys. Rev. Lett. 96, 197205 (2006). 

[3] Y. S. et ai, Phys. Rev. Lett. 91, 107001 (2003). 

[4] R. C. et al, Phys. Rev. Lett. 86, 1335 (2001). 

[5] S. Yan, D. A. Huse, and S. R. White (2011), 

arXiv:1011.6114. 
[6] Z. Y. Meng, T. C. Lang, S. Wessel, F. F. Assaad, and 

A. Muramatsu, Nature 464, 847 (2010). 
[7] X. Wen, Phys. Rev. B 44, 2664 (1991). 
[8] S. Sachdev, Phys. Rev. B 45, 12377 (1992). 
[9] T. Senthil and M. P. A. Fisher, Phys. Rev. B 62, 7850 

(2000). 

[10] O. Motrunich and T. Senthil, Phys. Rev. B 71, 125102 
(2005). 

[11] F. Wang, Phys. Rev. B 82, 024419 (2010). 

[12] Y.-M. Lu and Y. Ran, p. 1005.4229 (2010). 

[13] C. Xu, Phys. Rev. B 83, 024408 (2011). 

[14] L. Balents, M. P. A. Fisher, and S. M. Girvin, Phys. Rev. 

B 65, 224412 (2002). 
[15] D. Sheng and L. Balents, Phys. Rev. Lett. 94, 146805 

(2005). 

[16] S. V. Isakov, Y. B. Kim, and A. Paramekanti, Phys. Rev. 

Lett. 97, 207204 (2006). 
[17] S. V. Isakov, M. B. Hastings, and R. G. Melko, 

arXiv: 1102. 1721 (2011). 
[18] A. V. Chubukov, T. Senthil, and S. Sachdev, Phys. Rev. 

Lett. 72, 2089 (1994). 
[19] A. V. Chubukov, T. Senthil, and S. Sachdev, Nucl. Phys. 

B 426, 601 (1994). 
[20] S. V. Isakov, T. Senthil, and Y. B. Kim, Phys. Rev. B 

72, 174417 (2005). 
[21] A. Paramekanti, L. Balents, and M. P. A. Fisher, Phys. 

Rev. B 66, 054526 (2002). 
[22] O. I. Motrunich and M. P. A. Fisher, Phys. Rev. B 75, 

235116 (2007). 

[23] M. S. Block, R. V. Mishmash, R. K. Kaul, D. N. Sheng, 

O. I. Motrunich, and M. P. A. Fisher, Phys. Rev. Lett. 

106, 046402 (2011). 
[24] A. W. Sandvik, Phys. Rev. B 59, R14157 (1999). 
[25] H. Buchler, M. Hermele, S. Huber, M. P. Fisher, and 

P. Zoller, Phys. Rev. Lett. 95, 040402 (2005). 
[26] R. Melko and A. W. Sandvik, Phys. Rev. E 72, 026702 

(2005). 

[27] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. 

Fisher, Phys. Rev. B 40, 546 (1989). 
[28] M. Campostrini, M. Hasenbusch, A. Pelissetto, and 

E. Vicari, Phys. Rev. B 74, 144506 (2006). 
[29] T. Grover and T. Senthil, Phys. Rev. B 81, 205102 

(2010). 

[30] R. G. Melko, A. B. Kallin, and M. B. Hastings, Phys. 

Rev. B 82, 100409 (2010). 
[31] The q = (0, 0) peaks and their reflections only represent 

a finite density of plaquettes, not ordering 



